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摘 要 : 利用 最 新 的 GIMMS NDVI 3g + 数据 , 辅 以 MVC、 趋 势 分 析 、 变 异 系数 及 
9 然 文 化 遗产 保护 区 NDVI 时 空 分 布 特征 及 变化 趋势 ,并 对 其 驱动 因素 进行 分 析 。 研 
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首相 关 分 析 等 方法 ,探讨 了 


究 发 现 :@ 由 于 海拔 引起 的 水 热 条 件 的 变化 ,该 研究 区 NDVI 由 西北 向 东南 逐渐 增 大 ;@ 34 a 来 研究 区 NDVI 具有 
明显 增长 趋势 ,核心 区 NDVI 增 速 为 0.69% - (10a) ~! , 比 缓 冲 区 NDVI 增 速 高 出 0.41% * (10a) 5G) 研究 区 ND- 
VI 改善 面积 远大 于 退化 面积 ,总 体 上 呈 缓 慢 上 升 趋势 ,其 中 核心 区 NDVI 呈 改 善 和 退化 趋势 的 面积 分 别 占 47. 9296 
和 3.00% ,缓冲 区 中 有 65.70% 的 面积 基本 不 变 ; 由 34 a 间 研 究 区 NDVI 相对 较 低 波动 变化 和 中 等 波动 变化 共 占 全 


区 的 75.54% , 主 


分 布 在 核心 区 的 中 部 和 南部 以 及 缓冲 区 大 部 分 区 域 ;@) 相关 和 偏 相关 分 析 结 果 表 明 ,该 研究 区 


NDVI 受 温度 和 降水 共同 影响 ,与 降水 相 比 ,温度 是 影响 该 研究 区 NDVI 变化 的 主要 因素 。 


关键 词 : GIMMS NDVI 3g + ; 趋势 分 析 ; 


青藏 高 原 以 其 巨大 的 山体 效应 造就 了 亚洲 中 部 
腹地 的 干旱 核心 和 全 球 唯一 的 自然 地 域 格局 "”]。 
可 可 西里 世界 自然 文化 遗产 保护 区 位 于 青藏 高 原 北 
部 ,其 海拔 4 500 ~5 200 m, 低 温 少 雨 ,是 世界 著名 
的 寒 区 旱 区 及 高 寒 荒漠 过 渡 区 ,其 极端 脆弱 的 植被 
生态 环境 对 全 球 变化 响应 极为 敏感 。 

植被 是 连接 大 气 .水体 和 土壤 的 纽带 ,是 构成 陆 
地 生态 系统 的 主体 5 T 。 利 用 长 时 序 全 球 归 一 化 植 
被 指数 数据 集 的 波动 ,探视 全 球 化 气候 变化 规律 ,对 
评估 气候 变化 导致 的 生态 风险 ,研究 陆地 生态 系统 
的 演变 过 程 具有 重要 意义 。 利 用 GIMMS NDVI g 
和 CIMMS NDVI 3g 研究 植被 动态 响应 气候 变化 由 
来 已 久 。 例 如 ,Shen 等 利用 1982—2006 年 GIMMS 
NDVI g 数据 集 及 气候 资料 ,分 析 了 气候 变化 对 中 国 
草地 的 影响 “) 。Zhang 等 以 1982—2006 年 GIMMS 
NDVI g 数据 和 和 气象 数据 为 基础 ,对 气候 变化 下 的 青 
藏 高 原 植被 进行 了 分 析 ""”。 周 锡 饮 等 利用 GIMMS 
NDVI g 数据 对 蒙古 高 原 地 区 1981 一 2006 年 植被 覆 
盖 的 时 空 变化 进行 了 研究 ,并 从 气候 变化 和 人 类 活 
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动 的 角度 ,分 析 了 植被 覆盖 变化 的 原因 " 。 郭 鹏 
等 5 利用 1982—2006 年 的 GIMMS NDVI g 数据 研 
究 了 新 疆 的 植被 覆盖 动态 变化 情况 ,并 结合 气候 等 
相关 因素 综合 分 析 了 其 时 空 分 布 特征 。 与 老 版 本 的 
GIMMS NDVI g 数据 集 相 比 ,GIMMS NDVI 3g 数据 
集 在 检测 地 面 信息 方面 则 更 加 敏感 路 。 其 中 ,Du 
等 "4 使 用 全 球 库存 建 模 和 制图 研究 归 一 化 差异 植 
被 指数 (NDVI) 版 本 NDVI 3g 数据 集 , 分 析 了 青藏 高 
原 1982 年 以 来 14 个 渐 近 周期 的 NDVI 季节 分 量 趋 
势 及 其 与 气候 因子 的 关系 。 王 茜 等 利用 1982 一 
2012 年 的 GIMMS NDVI 3g 数据 ,研究 了 我 国 NDVI 
时 空 分 布 状况 及 空间 变化 率 , 同 时 选用 双 变 量 相关 
与 偏 相关 方法 ,分 析 了 我 国 不 同 植被 类 型 区 域 NDVI 
与 6 Pp T BTE TET. BEER AXE 
1982—2013 年 GIMMS NDVI 3g 数据 集 及 青海 50 个 
气象 观测 站 同期 气象 资料 ,利用 趋势 分 析 、 相 关 分 析 
方法 ,对 青海 省 植被 覆盖 时 空 变化 特征 及 其 驱动 因 
子 进行 了 分 析 "“。 昌 然 这 些 研究 利用 GIMMS ND- 
VI 数据 噪声 小 、 具 有 较 高 的 时 间 分 辨 率 的 特点 做 了 
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大 量 的 工作 ,但 受 其 数据 时 长 的 限制 ,未 能 对 相应 的 
研究 做 更 加 绩 密 的 分 析 。 

有 许多 学 者 关注 可 可 西里 生态 环境 的 变迁 ,但 
对 该 地 区 NDVI 变化 的 长 期 监测 、 空 间 分 布 `NDVI 
变化 与 气候 的 响应 及 人 类 活动 对 NDVI 变化 影响 等 
方面 的 研究 较 少 。 本 文采 用 1982 一 2015 年 最 新 的 
GIMMS NDVI 3g + 数据 集 ,结合 与 之 同期 的 气候 数 
据 ,对 可 可 西里 保护 区 长 时 序 NDVI 时 空 变化 及 气 
候 因子 进行 了 深入 的 研究 ,以 期 揭示 该 区 域 NDVI 
变化 的 时 空 分 布 和 变化 趋势 ,以 及 该 区 域 NDVI 变 
化 与 气候 响应 的 规律 ,可 为 世界 自然 遗产 地 的 保护 
及 管理 提供 参考 。 


1 研究 区 与 研究 数据 


1.1 研究 区 概况 

可 可 西里 蒙古 语意 为 "青色 的 山 染 ”, 位 于 青海 
省 玉树 自治 州治 多 县 、. 曲 麻 莱 县 境内 ,地 域 涵盖 可 可 
西里 国家 级 自然 保护 区 全 部 和 三 江 源 国家 级 自然 保 
护 区 部 分 。 绵 百 逾 500 km 的 昆仑 山 .可 可 西里 山 和 
乌 兰 乌 拉 山 勾勒 出 的 三 山 间 两 盆 , 核 心 区 面积 达 
3.73 x10? hm? ,缓冲 区 面积 2.29 x 10* hm? ,平均 海 
1&4 500 m 以 上 (图 1)。 可 可 西里 气候 干旱 寒冷 ， 
具有 典型 的 高 寒气 候 特 征 , 其 年 均 气 温 
-6.9 ~ -2.8 C ,年 均 降水 量 173 ~495 mm。 可 可 
西里 也 是 最 年 轻 的 高 原 冻 土 区 ,有 74 种 脊椎 动物 ， 
ABO E JEUEPEEAE . 藏 野 驴 RUSO , 狼 和 棕熊 的 栖 
息 地 ,其 中 生活 的 藏 羚羊 占 到 了 全 球 总 数 的 40% 。 
可 可 西里 地 处 青藏 高 原 高 寒 草 铅 高寒 荒 漠 的 过 渡 
区 ,主要 植被 类 型 为 高 寒 草原 和 高 寒 草 旬 ,分 布 有 较 
大 面积 的 高 山 冰 缘 植 被 。 由 于 严酷 的 气候 条 件 , 土 
Biévb pt ` 贫 将 ,植物 群落 结构 简单 ,生长 季 短 ,生产 力 
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Fig. 1 Geographical location and topography of Hoh Xil 


低下 ,植被 易 受 破坏 却 不 易 恢复 。 
1.2 研究 数据 

最 新 发 布 的 NetCDF 格式 的 GIMMS NDVI 数据 
是 NDVI 3g(1981 -07 -01—2012 -12 -31 ) 数 据 的 
加 长 版 ,来 源 于 ECOCAST 网 站 (https://ecocast. 
arc. nasa. gov/ data/pub/gimms/3g. v1/) ,时 序 延 长 至 
2015 年 ,时 长 达 35 a(1981 - 07 -01 一 2015 - 12 - 
31) ,空间 分 辨 率 为 0.083。x 0. 083" ,时 间 分 辩 率 为 
15 d? , 3 NDVI 数据 最 初 是 由 美国 航空 航天 局 
( National Aeronautics and Space Administration, 
NASA) 全球 监测 与 模型 研究 ( GIMMS ) 研究 组 采用 
国际 通用 的 最 大 值 合成 (maximum value composites, 
MVC) 法 合成 的 数据 集 '" ,投影 方式 为 Albers Coni- 
cal Equal Area 和 经 纬度 投影 两 种 ,其 原始 数据 源 于 
6 "i NOAA 卫星 (NOAA - 7,9,11,14,16,17) ,经 过 
辐射 校正 和 几何 粗 校 正 的 NOAA-AVHRR ,再 对 每 
H .每 轨 图 像 进行 几何 精 校正 \ 除 坏 线 \ 除 云 等 处 理 ， 
进而 进行 NDVI 计算 及 合成 “ ” 。 因 此 ,该 数据 相 
较 于 其 他 的 NDVI 数据 ,精度 更 高 .误差 更 小 ,适用 
于 植被 覆盖 变化 的 长 期 监测 。 本 文 所 使 用 的 数据 是 
经 过 IDL 程序 数据 格式 转换 .影像 裁剪 , 掩 膜 处 理 等 
预 处 理 过 程 ,将 原始 的 全 球 NetCDF 格式 数据 集 制 
作成 为 tf 格式 的 可 可 西里 地 区 数据 集 。 其 中 , 气 
象 数据 来 源 于 中 国 气象 科学 数据 共享 网 (http: / 
cdc. cma. gov. cn) 1982—2015 年 平均 气温 和 降水 
数据 。 


2 研究 方法 


2.1 NDVI 最 大 值 合成 
可 可 西里 世界 自然 文化 遗产 保护 区 处 于 青藏 高 
原 腹 地 的 高 海拔 地 区 ,NDVI 值 较 低 ,为 了 准确 反映 
该 地 区 NDVI 的 空间 分 布 及 变化 特征 ,采用 最 大 值 
合成 法 MVC 生成 月 NDVI 最 大 值 ,有 效 减 少 了 来 自 
云气 溶胶 、 云 阴影 、 视 角 以 及 太阳 高 度 角 的 影 
响 22-20 。 在 本 文 研 究 中 ,以 15 d 的 NDVI 数据 为 基 
础 ,采用 半月 合成 法 获取 月 度 最 大 值 NDVI, 公 式 为 : 
mNDVI, = Max NDVI 3g(a),NDVI3g(6)] (1) 


RP u 表示 月 序号 , 取 值 范围 为 1 ~ 12; mNDVIj x 
示人 第 j Hie K NDVI ff ;NDVI3g(a) ,NDVI3g(^) 
分 别 表示 第 j 月 上 半月 和 下 半月 的 NDVI 值 。 由 此 
获得 1 一 12 月 的 NDVI 年 内 数据 。 根 据 NDVI 年 内 
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数据 可 以 获得 植物 生长 的 周期 变化 。 
同样 采用 MVC 方法 ,在 月 度 最 大 值 合 成 的 基础 
上 对 NDVI 数据 进行 年 度 最 大 值 合成 , 其 公 
如 下 : 
aNDVI, = Max ( mNDVI, ,mNDVL ,--- ,mNDVI,, , 
(i 21982,1983 ,…2015 ) (2) 


从 而 得 到 NDVI 最 大 值 年 际 数据 ,该 数据 能 反 
映 NDVI 的 年 际 变化 规律 。 
2.2 ”NDVI 年 际 变化 趋势 

基于 像 元 的 NDVI 变化 趋势 ,通过 一 元 线性 回 
归 分 析 ,得 出 不 同时 间 NDVI 随时 间 变 化 的 趋势 , 即 
NDVI HERZL” ”® ,计算 公式 为 ; 


nx X(ixNDVL) - XiYNDVI, 


ope 7 


04 (3) 


nxë- (3i) 
式 中 :nm 表示 研究 的 累积 年 数 ;0 表示 年 最 大 ND- 
VI 在 n 年 的 变化 趋势 值 ;NDVI, 表示 第 i 年 研究 区 
年 最 大 NDVI 值 。 当 bope >0 时 ,表示 所 在 像 元 的 年 
最 大 NDVI 值 随时 间 旦 上 升 趋势 ; 当 Oae <0 时 , 表 
示 所 在 像 元 的 年 最 大 值 NDVI 随时 间 呈 下 降 趋势 ， 
M gas =0 时 ,表示 所 在 像 元 的 年 最 大 NDVI 未 随时 
间 变 化 。 
2.3 NDVI 变异 系数 

变异 系数 是 描述 随机 变量 离散 程度 的 统计 量 , 用 
以 分 析 NDVI 空间 格局 与 空间 分 异 规律 。 本 文采 
用 逐 像 元 计算 NDVI 变异 系数 的 方法 ,计算 公式 为 : 


(a-a)? 
a (4) 
式 中 :C, 表示 NDVI 的 变异 系数 ;n 表示 人 研究 的 累积 
年 数 ;x, 表示 第 i 年 年 最 大 NDVI; x 表示 1982— 
2015 年 年 最 大 NDVI 平 均值 。C, 值 越 小 ,表明 ND- 
VI 分 布 的 离散 程度 越 小 ;C, 值 越 大 ,表明 NDVI 分 
布 的 离散 程度 就 越 大 。 
2.4 NDVI 相关 与 偏 相 关系 数 

气候 变化 是 影响 该 研究 区 NDVI 变化 的 主要 因 
素 之 一 。 因 此 ,通过 逐 像 元 计算 年 最 大 NDVI 与 年 
均 气 温和 降雨 之 间 的 相关 系数 ,来 表征 气候 因子 与 
NDVI 变化 的 相关 程度 “1 。 其 计算 公式 为 : 


c-l 


X 


OX[ 696-5] 


Ry u n n 
VEC EG 


(5) 


ChinaXiv 合 作 期 刊 


Ur R 


RPR 表示 xy 两 变量 的 相关 系数 ;x; 与 y; 分 别 
表示 xy 两 变量 第 i 年 的 值 ;x* 和 了 分别 表示 两 变量 
n 年 的 平均 值 ;n 表示 研究 的 累积 年 数 。 在 简单 线 
性 相关 计算 的 基础 上 ,可 进一步 导出 偏 相关 系数 , 公 

cuc R, -RR, 

U XQ-RD X -R.) 
式 中 :R, .为 自 变 量 z 固定 后 因 变 量 * 与 自 变量 y 的 
偏 相 关系 数 。 


3 结果 与 分 析 


3.1 NDVI 空 间 分 布 特征 

参照 已 有 的 研究 ”0 ,对 34 a 可 可 西里 世界 自 
然 文化 遗产 保护 区 中 核心 区 和 缓冲 区 (以 下 分 别 简 
称 ”" 核心 区 ”和 "缓冲 区 ”) 年 最 大 NDVI 分 别 做 平均 
值 制图 ,得 到 研究 区 的 NDVI 分布 特 征 (图 2)。 

该 研究 区 的 西北 以 昆仑 山脉 为 界 ,地势 高 峻 , 东 
南 地 势 相对 平缓 ，。 如 图 2 所 示 :研究 区 NDVI 空 
间 分 布 的 差异 性 显著 ,青藏 公路 以 西 地 区 为 “可 可 
西里 藏 北 无 人 区 ?的 一 部 分 ,其 NDVI 值 在 -0.30 ~ 
0. 30 ,而 青藏 公路 以 东 地 势 平 缓 ,NDVI 值 相对 较 高 ， 
在 0.1~0.4。NDVI 由 西北 向 东南 随 海拔 的 降低 而 
逐渐 递增 ,该 结论 与 张 颖 等 ”的 研究 结果 一 致 ; 
NDVI 值 为 -0.30 的 像 元 为 水 域 , 以 深 色 表示 ,其 由 
西向 东 , 由 北向 南 ,在 核心 区 中 分 别 为 :可 可 西里 湖 、 
时 万 湖 .多 和 尔 改 湖 . 库 赛 湖 ,在 缓冲 区 中 分 别 为 : 西 金 
乌 兰 湖 . 乌 兰 乌拉 湖 。 
3.2 NDVI 年 际 变化 规律 

分 别 对 核心 区 和 缓冲 区 1982—2015 年 最 大 
NDVI 均值 进行 分 析 ,获得 年 最 大 NDVI 均值 分 布 结 
果 ( 图 3)。 
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图 例 F 
| 口 核心 区 国 冰 川 和 永久 积 雪 C 
缓冲 区 加 沼泽 湿地 DAREM 0 50km 
mk mji 咽 高 寒 草 旬 
90* E 9PE 922 E 93*E 94^ E 
图 2 研究 区 34 a 年 最 大 NDVI 平 均值 空间 分 布 
Fig.2 Spatial distribution of the average of annual maximum 
NDVI in 34 years 
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一 -核心 区 NDVI 均 值 ”一 核心 区 线性 趋势 
024p ”一 缓冲 区 NDVI 均 值 ”一 缓冲 区 线性 趋势 
0.22 | s 

im N A M 
g * 
E 0.20 
^ 
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0.16 L 


1980 1985 1990 1995 2000 2005 2010 2015 
年 份 


图 3 研究 区 34 a 年 最 大 NDVI 均值 变化 趋势 
Fig.3 Change trend of the average of annual maximum 


NDVI in the study area in 34 years 


由 图 3 可 知 :核心 区 年 最 大 NDVI 均值 比 缓冲 
区 低 0.037 0; 年 最 大 NDVI 从 1982 一 2015 年 整体 
呈 上 升 趋势 ,其 中 核心 区 年 最 大 NDVI 均值 在 34 a 


间 的 增 速 为 6.930 5 x 10 * a^! ,而 缓冲 区 年 最 大 
NDVI 均值 的 增 速 为 2.782 6 x 10 *a^! Bp t DX 4E 
最 大 NDVI 的 增长 速率 约 是 缓冲 区 的 2.49 fi. 34 a 
间 核 心 区 与 缓冲 区 的 NDVI 均值 波动 基本 一 致 ,最 
小 值 均 出 现在 1995 年 ,分别 为 0.156 5 和 0. 191 8; 
最 大 值 均 出 现在 2010 年 ,分 别 为 0.203 1 和 0.235 7。 
其 中 ,1982 一 2000 年 整体 呈 上 升 趋势 ,2000 一 2007 
年 整体 呈 下 降 趋 势 ,2007 一 2010 年 快速 上 升 至 最 大 
值 ,之 后 持续 下 降 至 2015 年 。 
3.3 NDVI 年 际 变 化 分 析 
根据 公式 (3 ) ,通过 一 元 线性 回归 对 年 最 大 

NDVI 数 据 进行 逐 像 元 分 析 , 获得 逐 像 元 回归 曲线 
斜率 9, 图 像 ( 图 4)。 参 照 已 有 的 研究 ,根据 

bw 将 研究 区 中 的 NDVI 变化 分 为 5 个 等 级 , 即 严重 
退化 ,轻微 退化 、 基 本 不 变 、 轻 微 改 善 和 明显 改善 
( 表 1)。 


表 1 34 a NDVI 年 际 变 化 趋势 统计 
Tab.1 Interannual variation trend of NDVI in 34 years 


变化 情况 bwe 变 化 范围 核心 区 比例 /% 核心 区 面积 /km 缓冲 区 比例 /% 缓冲 区 面积 /km 
严重 退化 buoe< -0.002 0 0. 64 236.43 0.79 182.07 
轻微 退化 -0.0020 «0,4, S -0.000 7 2.36 874.68 9.23 2 124.01 
基本 不 变 -0.0007 ah «0.000 7 49.20 18 204. 48 65.70 15 110.81 
轻微 改善 0.000 7 < 64, 0. 002 0 44. 09 16 313.08 22.69 5 218.00 
明显 改善 baone >0.002 0 3.83 1 418.54 1.58 364.11 

总 计 100 37 000 100 23 000 


吕 核心 区 — A 
- 缓冲 区 — 国 基 本 不 变 0 50km 


m 严重 退化 国 轻 微 改善 ”国明 显 改善 一 一 
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图 4 研究 区 1982—2015 年 NDVI 年 际 斜率 变化 
Fig.4 Variation of interannual slope of NDVI in the study 
area during the period of 1982 — 2015 


分 析 可 得 :(1)NDVI 呈 增长 趋势 的 有 23 314. 73 
km? , 占 研究 区 面积 的 36. 10% ; 呈 基 本 不 变 趋势 的 
有 33 315.30 km , 占 研究 区 面积 的 37.45% ;时 减少 
趋势 的 有 3 417.18 km , 占 研究 区 面积 的 6.51% 。 

(2) 研 究 区 中 NDVI 年 际 变化 呈 增 长 趋势 的 面 
积 约 是 呈 减 少 区 域 面 积 的 6. 82 倍 ; NDVI 年 际 变化 


呈 增 长 的 区 域 主要 分 布 在 核心 区 的 中 部 和 北部 ;组 
冲 区 的 大 部 分 区 域 其 年 际 变化 趋势 基本 不 变 
3.4 NDVI 变异 系数 分 析 

根据 公式 (4) , 逐 像 元 计算 1982 一 2015 年 核心 
区 及 缓冲 区 NDVI 的 变异 系数 C, ,并 根据 C, 的 大 小 
将 其 分 为 5 PER : 低 波 动 变化 .相对 较 低 波动 
变化 .中 等 波动 变化 .相对 较 高 波动 变化 和 高 波动 变 
化 ( 表 2) ,核心 区 和 缓冲 区 NDVI 变异 系数 如 图 5。 

由 图 5 可 以 看 出 ,变异 系数 反映 的 NDVI 波动 
规律 是 :相对 较 低 波动 占 比 最 高 ,中 等 波动 次 之 , 低 
波动 主要 为 水 域 ;高 波动 相对 较 高 波动 及 中 等 波动 
的 区 域 主要 分 布 在 核心 区 的 北部 ;相对 较 低 波动 的 
区 域 主要 分 布 在 研究 区 的 东南 部 ; 低 波动 变化 的 区 
域 所 占 面积 最 小 ,其 位 置 分 别 对 应 着 湖泊 和 雪山 ,地 
形 地 貌 对 NDVI 的 离散 程度 具有 一 定 的 影响 。 
3.5 NDVI 与 气候 因子 相关 性 分 析 

在 全 球 变 暖 的 背景 下 5 ,研究 区 1982 一 2015 
年 的 增 温 趋 势 显著 (图 6a 和 图 6b) ,与 此 同时 ,研究 
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图 例 


国 相 对 较 低 波动 变化 。 国 高 波动 变化 
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图 5 


1982—2015 年 研究 区 多 


FE 最 大 NDVI 变异 系数 


Fig.5 Variable coefficient of annual maximum NDVI in the 


study area during the period of 1982 — 2015 


R2 研究 区 34 a 年 最 大 NDVI 的 变异 系数 
Tab.2 Variable coefficient of annual maximum 


NDVI in the study area in 34 years 


ChinaXiv 合 作 期 刊 


区 的 降水 量 整 体 也 呈 增 加 趋势 (图 6c、 图 6d)。 由 图 
6 可 得 ,核心 区 年 平均 气温 的 上 升 速率 为 0.056 © - 
a , 比 缓冲 区 年 均 气 温 的 上 升 速率 低 0.01 C «75; 
而 核心 区 年 降水 量 的 增加 速率 为 3.67 mm a”, 比 
缓冲 区 年 降水 量 的 增加 速率 高 出 0.25 mm * a^, 

依据 公式 (5 ) , 逐 像 元 计算 1982—2015 年 年 最 
X NDVI 与 年 均 气 温和 降水 的 相关 系数 ,获得 相关 
系数 图 (图 7 和 图 8) ,分 析 可 得 : 

(1) 年 最 大 NDVI 与 年 均 气 温和 降水 的 相关 系 
数 分 别 为 0.248 和 0.191, 其 中 ,核心 区 NDVI 与 其 
气温 和 降水 的 相关 系数 分 别 为 0.303 和 0. 230 ,组 
冲 区 NDVI 与 其 气温 和 降水 相关 系数 分 别 为 0. 158 
fill 0. 126. 

(2) 年 最 大 NDVI 与 气温 和 降水 成 正 相 关 的 区 
域 主要 分 布 在 核心 区 的 中 部 和 北部 以 及 缓冲 区 南部 
地 区 ,其 核心 区 和 缓冲 区 中 NDVI 与 年 均 气 温和 降 
水 量 的 相关 性 统计 见 表 3。 

进一步 对 年 最 大 NDVI 与 年 均 气 温和 降水 量 进 
行 偏 相关 分 析 , 获得 其 偏 相关 系数 图 (图 9 和 图 


10) ,分 析 可 得 : 


(1) 年 最 大 NDVI 与 气温 和 降水 的 平均 偏 相关 


本 ALES 占 核心 区 ” 占 缓 冲 区 
变异 程度 Cv 的 变化 范围 面积 比例 /% 面积 比例 /% 
低 波动 变化 Cv « 0. 05 8.31 5.63 
相对 较 低 波动 变化 0.05 < Cv<0. 10 36. 10 62. 54 
中 等 波动 变化 0.10 < Cv<0. 15 32.43 20. 00 
相对 较 高 波动 变化 0.15 < Cv<0.20 13. 10 7.04 
高 波动 变化 Cv >0. 20 10. 06 4. 79 
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-3.5 r (a) 核心 区 0.21 
e 
s z 
w 
E 三 
Hu 
1980 1985 1990 1995 2000 2005 2010 2015 
一 -年 均 降水 
450 [ (QO 核心 区 - 0.21 
400 f 4 . AA. 1030 
E E "E T m 0.19 
ug 300 . A " j 
pi Mi MS 018 Z 
WE 250] W Vy y 2 
A 200} a / 2. i ii 0.17 
150[ M 0.16 
100 nr 


一 -年 最 大 NDVI 均 值 
-2.5 [ (b) 缓冲 区 1024 
-3.0 ; M i 
-3.5 T. A n | ues 
p -40 x Ji ( "NI yi 10.22 
如 -5.0 上， f | 人 021 2 
Sob ` n" 
-60 上 | ino 
-6.5 10.19 
-7.0 
E EEN PEE E PERLER E: 
1980 1985 1990 1995 2000 2005 2010 2015 
一 -年 最 大 NDVI 均 值 
500 | (ay 缓冲 区 40.24 
ad i À , 10.23 
g 40- , j FA A 
E N^ "IR m ^y 10.22 
* NAR VINA ozn & 
& 300 | A^ /RY V ) 1021 B 
ri om ~ VV 10.20 
200 | 
i 10.19 
100 0 
1980 1985 1990 1995 2000 2005 2010 2015 


1980 1985 1990 1995 2000 2005 2010 2015 


图 6 


1982—2015 年 研究 区 年 NDVI 均值 与 气候 因子 的 相关 性 


Fig.6 The correlation between NDVI and climate factors in the study area during the period of 1982 -2015 
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陈 登 魁 等 :1982 一 2015 年 可 可 西里 NDVI 变化 特征 及 其 气候 响应 


表 3 研究 区 NDVI 与 气候 因子 的 相关 性 统计 
Tab.3 Correlation between NDVI and climatic 
factors in the study area 

NDVI 与 年 均 气 温 NDVI 与 年 均 降水 
相关 程度 核心 区 缓冲 区 核心 区 ”缓冲 区 


比例 /% ”比例 /% 比例 /% ”比例 /% 
非常 显著 负 相 关 1.12 3.96 0. 64 20.32 
显著 负 相 关 7.51 18.21 13. 42 1.06 
不 相关 8.79 12. 12 6.23 15.02 
显著 正 相关 35. 14 39.31 47.60 47.50 
非常 显著 正 相 关 — 47.44 26.40 32. 11 16. 10 


À 


Z) 口 缓冲 区 BRAVE 
e| 国 非 常 显著 负 相关 国 显 著 正 相关 0 50km 


国 显 著 负 相关 有 非常 显著 正 相 关 = 
90°E 91°E 92°E 93°E 94°E 
图 7 研究 区 NDVI 与 气温 的 相关 性 


Fig.7 The correlation between NDVI and temperature in 


the study area 
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图 8 研究 区 NDVI 与 降水 量 的 相关 性 


Fig.8 The correlation between NDVI and precipitation in 


94^ N 


the study area 


系数 分 别 为 0.420 和 0.311 ,其 中 核心 区 NDVI 5E 
气温 和 降水 的 平均 偏 相关 系数 分 别 为 0. 523 和 
0.374 ,缓冲 区 NDVI 与 其 气温 和 降水 的 平均 偏 相 关 
系数 分 别 为 0.257 和 0.207。 

(2) 年 最 大 NDVI 与 气温 和 降水 呈正 相关 区 域 
分 别 占 研究 区 的 79.20% 和 74. 93% ,其 中 核心 区 和 
缓冲 区 中 NDVI 与 气温 成 正 相 关 的 区 域 分 别 占 
52. 64% 和 26.57% ,而 与 降水 成 正 相 关 的 区 域 分 别 
占 49.06% 和 25.87% 。 


36° N 


图 例 


口 核心 区 
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9 研究 区 NDVI 与 气温 的 偏 相关 性 
Fig.9 The partial correlation between NDVI and 


temperature in the study area 
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到 10 研究 区 NDVI 与 降水 的 偏 相关 性 
Fig.10 The partial correlation between NDVI and 
precipitation in the study area 
4 结论 


(1) NDVI 均值 年 际 分 析 表 明 :研究 区 NDVI 的 
空间 分 布 特征 与 高 海拔 影响 下 的 水 热 特 征 高 度 一 
致 , 即 从 东南 向 西北 海拔 逐渐 升 高 ,气温 逐渐 降低 而 
降水 逐渐 减少 ,其 植被 类 型 则 体现 出 由 “高 寒 草 甸 ” 
DM EE SUN MEUSE S MEN E SUE 

(2) NDVI 年 际 变化 规律 :1982 一 2015 年 研究 
区 NDVI 在 34 a 间 整体 呈 增 长 趋势 ,但 有 波动 。 其 
间 ,在 1995 年 出 现 最 小 值 ,在 2010 年 出 现 最 大 值 ， 
其 中 核心 区 NDVI 增 速 为 0.69% - (10a) ^! , EZ up 
[X NDVI 增 速 高 出 0.41% . (10a)-。 

(3) 趋势 分 析 表 明 : 研 究 区 NDVI 改善 面积 远 
大 于 退化 面积 ,总 体 上 呈 缓 慢 改 善 趋势 ,其 中 核心 区 
NDVI 呈 改 善 和 退化 趋势 的 面积 分 别 占 47. 9296 和 
3.0096 ,缓冲 区 中 有 65.70% 的 面积 基本 不 变 。 

(4) 变异 系数 分 析 表 明 :34 a 间 研 究 区 NDVI 
相对 较 低 波动 变化 和 中 等 波动 变化 共 占 全 区 的 
75.54% ,主要 分 布 在 核心 区 的 中 部 和 南部 以 及 缓冲 
区 大 部 分 区 域 。 
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(5) 气候 相关 性 分 析 表 明 :NDVI 与 年 均 气 温 
的 相关 性 相 较 于 年 均 降水 ,其 相关 程度 更 高 , NDVI 
与 年 均 气 温 的 平均 相关 系数 比 NDVI 与 降水 的 平均 
相关 系数 高 0.057。 偏 相关 分 析 同 样 表 明 ,NDVI 与 
年 均 气 温 的 相关 程度 高 于 降水 。 可 见 ,温度 是 影响 
研究 区 NDVI 变化 的 主要 因素 。 
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Abstract: In this study ,the latest data of GIMMS NDVI 3g plus, MVC ,trend analysis , coefficient of variation and 
partial correlation analysis were used to study the spatial and temporal distribution and change trend of NDVI in the 
Hoh Xil Natural and Cultural Heritage Reserve in Qinghai Province during the period of 1982 — 2015 and to analyze 
the driving factors. The results showed that; (I) The NDVI in the study area increased gradually from northwest to 
southeast due to the change of water and heat conditions caused by altitude; (2) During the period of 1982 —2015, 
the NDVI in the study area was in an obvious increase trend ,the increase rate of NDVI in the core was 0. 69% ， 

(10a) ^! ,and it was 0. 4196 - (10a) '' higher than that in the buffer area; (3) The area of the regions with NDVI 
improvement in the study area was much larger than that with NDVI degeneration , and the NDVI was holistically in 
a slow increase trend. The areas of the regions with NDVI improvement and NDVI degeneration in the core accoun- 
ted for 47.92% and 3. 0096 respectively ,65. 7096 area in the buffer area was basically maintained; (4) The area of 
the regions with relatively low and moderate NDVI fluctuation in the 34 years accounted for 75. 54% of the total, 
and these regions were mainly distributed in the central and southern parts of the core and most areas of the buffer 
zone; (5) The results of correlation and partial correlation analysis showed that the NDVI in the study area was af- 
fected by temperature and precipitation , especially by temperature. 

Key words; | GIMMS NDVI 3g Plus; trend analysis; partial correlation analysis; Hoh Xil World Natural and Cul- 


tural Heritage Reserve; 


